Synergy between surfactants’ stiffness and concentration on their self-assembly into reverse micelles as water droplet carriers in nonpolar solvents

A study of the self-assembly process into reverse micelles (RMs) of linear surfactants and monomeric aqueous solutes dissolved in nonpolar solvents, varying the concentration (cs) and the persistence length (Lp) of the surfactants is presented here. The influence of cs and Lp on the structural and dynamic properties of the aggregates is investigated through mesoscopic simulations carried out with the dissipative particle dynamics method. All simulations are performed at a fixed water/surfactant molecular ratio of 2:1, varying the surfactant concentration from c = 6 wt% up to c = 12 wt%, for increasing surfactants’ rigidity from Lp = 0.73 nm up to Lp = 44.99 nm. It is found that there exists a collaborative interplay between cs and Lp that enhances the number of RMs assembled and their diffusion as carriers of water droplets. These results should be useful as guidelines to understand and improve processes where the RMs are implemented to carry aqueous solutes in nonpolar solvents.


Introduction
Reverse micelles (RMs) are self-organized agglomerations formed by amphiphilic surfactants when these are dissolved in nonpolar or organic solvents.Although there is evidence indicating that the RMs formation is influenced by the addition of water [1,2] or polar solutes [3,4], the self-assembly of surfactants into stable RMs in nonpolar solvents also occurs in the absence of polar solutes [5,6].RMs have a wide range of technological applications, such as in the synthesis of nanomaterials [7][8][9], catalytic processes [10,11], and drug delivery systems [12,13].They are important also in different industrial processes, such as in the food industry as protein extractors [14,15], and as thickeners of supercritical carbon dioxide (scCO 2 ) for the enhanced oil recovery and fracking industries [16][17][18][19].
Within the scope of numerical simulations, the atomistic study of reverse micellar systems has provided useful insights about the mechanisms behind the self-assembly process at molecular-and atomistic-levels [20][21][22][23][24][25].Earlier molecular dynamics (MD) simulations carried out by Brown and Clarke [21], provided information about the structural properties of a single insolated RM as a water droplet carrier in a nonpolar solvent.On the other hand, studies of the self-assembly process forming RMs in specific nonpolar solvents arose a decade after.For instance, Tobias and Klein [26] were the first to investigate the local properties of RMs formed by a mixture of water/calcium carbonate/calcium sulfonate, in the nonpolar carbon tetrachloride and octane solvents.Salaniwal and collaborators [22] were the first to study the selfassembly process of water/surfactant [(C 7 F 15 )(C 7 H 15 )CHSO 4 -Na + )] aggregates in scCO 2 .Allen et al. [23] studied the local properties of RMs formed by water/poly(oxyethylene) surfactants in decane as solvent.
The importance of coarse-grained numerical methods relies on their capacity to simulate systems at larger length and time scales in comparison with atomistic-level techniques, such as MD.Dissipative particle dynamics (DPD) is a coarse-grained simulation method developed to study the hydrodynamic properties of complex fluids [27,28].Several DPD works have delved into the study of the mechanisms that lead to the self-assembly of different classes of surfactants into direct [29,30] or reverse micelles [31][32][33].Recently, Lavagnini et al. [34] provide a detailed report about the parametrization of the chemical structure of various types of ester-, amide-, and sugar-based surfactants into DPD-parametrized structures.They modeled the self-assembly of 27 different surfactants in aqueous solvents and found that the critical micelle concentration predicted by their simulation agreed with the available experimental data [34].The work carried out by Lin and co-workers [35] directly implemented a RM-based system for the extraction of proteins, using the DPD method.They proved that results obtained for the recovery efficiency of the protein were validated with experimental data found in the literature [35].
Despite all the work there is around RMs, it is still necessary to invest efforts in the study of these systems to improve the understanding of both the self-association process in the formation of aggregates and the technological/industrial implementation of RMs.This work presents the study of structural and dynamic properties of the self-assembly process of linear surfactants into RMs in nonpolar solvents [33].It has been reported that the self-assembly of linear surfactants into RMs is modulated by their persistence length (L p ), for fixed molecular water/surfactant ratio of 1:1 and at constant surfactant's concentration (c s ) [33].Here, cooperative effects are reported on the structural and dynamic properties of the self-assembly process forming water/surfactant aggregates in nonpolar solvents.Our results show that all RMs have roughly spherical shape regardless of the concentration or the surfactant's rigidity.However, the structure of the surfactant's tail groups that envelopes the water/surfactant headgroups' nuclei depends on the rigidity degree of the surfactants.It is also found that the number of RMs and their diffusion depend on the competition/cooperation between the surfactants' concentration and their persistence length.

(a) Dissipative particle dynamics method
All simulations were carried out using the coarse-grained DPD simulation method [27,36].The DPD model is similar to the standard MD simulation algorithm, in the sense that the equation of motion of all the particles making up the system is solved simultaneously to obtain their position and momenta [37].The DPD force-field consists of three pairwise, additive central forces, namely the conservative (F C ij ), dissipative (F D ij ) and the random forces (F R ij ) [36].The latter two forces are given by where r ij and v ij are the relative position and velocity vectors between the i-th and j-th particles, respectively; and rij ¼ r ij =jr ij j is the unit position vector.The symbols ξ ij = ξ ji are random numbers with Gaussian properties.The weight functions o make of these short-range forces that drop to zero when the relative distance between the i-th and jth particles (r ij ) is larger than the cutoff radius (r * c ).The parameters γ and σ are the maximum strength of the dissipative and random forces, respectively, which are coupled through the fluctuation-dissipation theorem.The balance between these forces yields a built-in thermostat, such that k B T = σ 2 /2γ [36].Here, k B is Boltzmann's constant and T is the absolute temperature.
The remaining of the three fundamental DPD forces is the conservative force (F C ij ), defined as a soft, linearly decaying repulsive force that also drops to zero when r ij � r * c : In Eq (3), a ij is the maximum force intensity between the i-th and j-th particles.Table 1 shows the values of the parameters a ij used in this work, which have been successfully tested in previous reports.In particular, with the parameters in Table 1, faster aggregation was achieved of both direct (using Gemini surfactants) [38] and reverse micelles (using linear surfactants) [33] in aqueous and nonpolar solvents, respectively.Further details for the calculation of a ij parameters can be found elsewhere [39].
To model chain-like, (H)ead-(T)ail linear structures, such as the HT5 surfactant, with one DPD head bead and 5 beads in the tail used in this work, two harmonic, conservative forces are added to the model.One of them models the bonding forces between neighboring surfactant beads.The Kremer-Grest bead-spring model [40] is used to join two consecutive beads along the surfactant, using a harmonic oscillator: where k S and r 0 are the spring constant and its equilibrium distance, respectively.The second harmonic force is introduced to model the rigidity of each surfactant chain, which is an angular spring based on a three-body force, given as follows: where θ ijk is the angle formed by three consecutive beads in the surfactant.In Eq (5), k A and θ 0 are the angular spring constant and the equilibrium angle, respectively.

(b) Simulation details
All quantities used in the simulations are expressed in reduced DPD units, which are listed in Table 2.
All simulations were performed under canonical ensemble conditions, where the number density and the global temperature of the systems are always fixed.All systems studied here were placed in a cubic simulation of L * x ¼ L * y ¼ L * z ¼ 17:1r * c ; using periodic boundary conditions in all the faces of the box.The total number of DPD particles in all systems was fixed at N = 1.5 × 10 4 , to keep a constant reduced number density of ρ * = 3.The reduced temperature was always kept equal to one [36].All simulations reported here were carried out in two stages.In the first, systems were run for over 30 blocks of 10 4 time steps (Δt) each, to allow the system to reach thermal equilibrium.In the production stage, the systems were run for an additional 20 blocks of 10 4 Δt each, where most of the properties of interest are averaged.Considering the two stages, the total simulation time of each system was *45.18 ns.4)) and F A ij (Eq ( 5)), respectively.The equilibrium distance and angle are r 0 and θ 0 ; see F S ij (Eq (4)) and F A ij (Eq (5)), respectively.In (c) and (d), one finds simplified representations of the aggregates formed by flexible and rigid surfactants, respectively.https://doi.org/10.1371/journal.pone.0294913.g001Table 2. Reduced DPD units and their conversion factors.The SI conversion factors correspond to the scenario when the coarse-grained degree is N m = 3, which means that each DPD particle represents a group of three water molecules at room temperature [39,41].The asterisks refer to reduced DPD units.The parameters of the Kremer-Grest bead-spring force (Eq (4)) of each HT5 surfactant were chosen as k S ¼ 100k B T * =r * 2 c and r 0 ¼ 0:7r * c [33].These parameters have been successfully used to model different types of linear-chain structures [33,42], demonstrating that are appropriate to prevent bond crossing [43].The interplay between surfactant concentration and persistence length on the water/surfactant self-assembly into RMs, was studied for four surfactant concentrations (c s ) and five values of the persistence length (L p ).Here, the value of the angle of equilibrium was set in θ 0 = 180˚to model a linear structure, while the persistence length is related to the angular spring constant k A (see Eq (5)) as follows [44]:

DPD units SI conversion factors
In Eq ( 6), b 0 is the average bond length, such that b is the Langevin function [44].Table 3 shows the different surfactant concentrations and persistence lengths modeled in this work.

Results
To understand the role played by the concentration of linear surfactants and their rigidity on their self-assembly into RMs, a set of simulations of water/HT5 surfactant mixtures at fixed 2:1 molecular ratio was carried out.The surfactant concentration was changed proportionally to the total weight percentage of the system, considering four different surfactant concentrations, c s .The rigidity of the HT5 surfactant was given by its persistence length, L p (Table 3 lists the values of c s and L p studied in this work).Thus, the cooperative effects of c s and L p on the structural and dynamic properties of the water/HT5 aggregates formed in nonpolar solvents, were studied from a mesoscopic point of view.
Let us start with the structural analysis of the density profiles of the systems shown in As expected, the self-assembly process yields the formation of RMs, and the water beads dissolved in the nonpolar solvent are encapsulated in the RMs.It is to be noted that in all cases, all the water beads dissolved in the system are captured inside the micelles' nuclei, as seen in the snapshots in each panel of Fig 2 .The blue and red particles correspond to the water beads and surfactants' headgroups, respectively.This encapsulation can also be observed in the density profiles, since the water beads' profiles (blue lines in   Recently, it has been demonstrated that the number of RMs formed by linear surfactants models, such as the one used in this work, can be measured from their density profiles by counting the number of maxima of water beads or surfactants' heads (components of the RMs nuclei) [33].However, in the surfactant concentration range studied here (from c s = 6 wt% up to c s = 12 wt%), it becomes more difficult to identify how many aggregates each maximum corresponds to.For instance, in the case of the most rigid surfactants (L p = 44.99nm) at c s = 6 wt %, two RMs are formed due to the water/HT5 aggregation (see snapshot in  Different simulation studies have demonstrated that the number of RMs formed by linear surfactants in nonpolar solvents depends on the system's conditions.Among those are the water/surfactant molecular ratio and the type of solvent in which the micellization process occurs [2].The number of RMs depends also on some physical properties of the surfactants, such as their stiffness and tail group's length [33].Our results show that by increasing c s at fixed water/HT5 molecular ratio of 2:1, higher number of RMs are found, which encapsulate all water beads dissolved in the nonpolar medium, regardless of the stiffness of the HT5.On the other hand, fewer RMs are formed when increasing L p for a fixed surfactant concentration in most cases, as seen in Fig 3, where the average number of RMs (N ave RMs ) is plotted as a function of L p for the four different concentrations.
As Fig 2 shows, the HT5 surfactant model used in this work exhibits perfect efficiency as water-droplet container.In all the systems studied, either at different surfactant concentrations or for the different surfactant's persistence length, all the water beads dissolved were encapsulated by the RMs, see also S1 Fig in the S1 File.Furthermore, it is found that the self-assembly process leading to the RMs formation is stable as a function of time, since the number of RMs formed did not change throughout the course of the simulations.This is in contrast to the case of Gemini surfactants, whose self-association process leads to the formation of unstable micelles, with dynamic aggregation/segregation process for different types of architectures such as P-shaped [45] or X-shaped [38].The aggregation mechanism leading to the formation of RMs reported here for the HT5 surfactant is stable enough that if two RMs fuse into a single RM, no water beads inside the nuclei are released during the aggregation process.
To demonstrate the stability of the micellization process, the time evolution of the number of RMs (N RMs ) formed by the most flexible (Fig 4 These results also reveal that there exists synergy between L p and c s , leading to enhanced self-assembly between the water beads and the HT5 surfactant.The number of aggregates formed during the micellization process depends on L p and c s ; thus, the more rigid the surfactants are, the lower the number of aggregates is.On the other hand, the number of aggregates tends to grow with surfactant concentration, regardless of their rigidity. Having shown that the HT5 surfactant model used here leads to the formation of stable RMs, now the analysis of the structural features of the water/HT5 aggregates is presented.This is done through the calculation of the averaged radial distribution functions (RDFs), g(r), for specific pairs of components of the systems.The RDFs allow one to understand the effects produced by L p and c s on the formation of RMs and provide information about the local structure of the system.(W-H) RDFs (r = 0.56 in Fig 5(A) and 5(B)).These structural characteristics are typically obtained when modeling RMs with linear HT surfactant models, as reported by Mayoral et al [2].who study the water/HT4 self-assembly into RMs in nonpolar solvents.
The influence of surfactant concentration on the RDFs between the water beads and the surfactant head-and tail-groups, Fig 5(A)-5(D), respectively, is analyzed at short and large relative distances.At short distances, the spatial correlation between water beads and HT5 surfactant chains shows that the tail groups surround the RMs nucleus made up of the water beads and head groups.On the other hand, both (W-H) and (W-T) RDFs show weaker correlation as c s increases, because the average intensity of g(r) decreases with the increasing concentration of water beads and surfactant chains in the system.Additionally, the RDFs in Fig 5 reveal that the number of aggregates formed in the system affects directly the spatial correlation between water beads and HT5 at relatively large distances.That is because at large distances, the water beads in the RMs nucleus only correlate with the surfactants' head-or tail-groups belonging to other aggregates.Therefore, the (W-H) and (W-T) correlations tend to have different decaying rates at large relative r values, depending on the number of RMs.It is typically found that the (W-H) and (W-T) RDFs tend to decay at relatively short distances if the number of aggregates is higher; see Fig 5(A) and 5(C).It is remarked that these effects occur because in this range of concentrations (6 wt% � c s � 12 wt%), the water beads/HT5 aggregates formed are always RMs-like structures, with no free water beads found in the system.
The influence of the surfactant's persistence length on the structure of the micelle nuclei is minimal at short distances, as seen in the (W-H) RDFs in Fig 5(A) and 5(B).Focusing on large relative distances, say at r > 5, the (W-H) correlations only occur between different micelle's nuclei.Therefore, the (W-H) g(r)'s tend to unit at large relative distances because in all cases exists more micelles surrounding the (W-H) beads.The consequences of increasing of L p on the (W-T) spatial correlations are shown in Fig 5(C) and 5(D), corresponding to the most flexible and most rigid surfactants, respectively.Firstly, it is noted that the intensities of the (W-T) g(r)'s decrease as the surfactants gain rigidity, since for stiffer surfactants the aggregation number grows; see also Fig 4. On the other hand, it is found that the (W-T) spatial correlation prevails over further distances for the stiffer surfactants because the tail groups are more stretched around the RMs nuclei.
Next, the structure surrounding the RMs nuclei was studied through the analysis of the spatial correlation between the surfactants' head-and tail-groups.To gain further understanding of the HT5 structure in the RMs, the mean end-to-end distance R e h i of HT5 surfactant is calculated, as follows: where r l n and r f n are the position vectors of the last (l) and first (f) particles of the n-th surfactant chain, respectively.The vertical bars, | |, refer to the magnitude of the vector r l n À r f n , and the angular brackets represent the average over time.Table 4 shows the mean end-to-end distances of the HT5 surfactants.
The results listed in Table 4 show that R e h i increase as the surfactant persistence length does, tending to their fully extended length (R FEL ) as the surfactant stiffness increases.To define the fully extended length of the surfactant, the average bond distance of adjacent surfactant's beads was used.In this case, such bond distance is found to be b 0 h i ¼ 0:76r * c , indicated by the first maximum in Fig 6(C).Thus, the fully extended length of the surfactant is defined as: with N b as the total number of bonds in the surfactant.
The results obtained for R e h i as a function of persistence length are shown in Fig 7, where the average of R e h i is reported over the four concentration values, since R e h i does not show significative differences between the c s values (see Table 4).In Fig 7, our predictions are compared with those reported by Herna ´ndez et al. [33] for larger, HT9 surfactant chains.The  results of Firetto et al. [46], who study the effects of the stiffness on the micellization of the H4T4 surfactant model, are included as well.
Additionally, the dynamics of the RMs are studied by means of their self-diffusion coefficient, D, as a function of the surfactants' persistence length, for each of the surfactant concentrations.The self-diffusion coefficient D is calculated from the mean-square displacements (MSD) of the center of mass of the RMs, during the last 10 simulations blocks, which is equivalent to 1 × 10 5 Δt � 9 ns, through Einstein's relation [47] MSD The self-diffusion coefficients of the RMs are shown in Fig 8 .The values of D reported in this work are comparable to those of RMs used to encapsulate water molecules in experiments [48,49].It is also found that the combined influence of the surfactants' concentration and persistence length can improve the mobility of water beads/HT5 aggregates, in the regime of lower concentration, for flexible surfactants.Our results show that, in this latter regime, the diffusion capacity of the RMs is low because a higher number of aggregates are formed; therefore, the RMs find it more difficult to move freely.However, when the aggregates are assembled by semiflexible HT5 surfactants, their diffusion is enhanced, although the concentration dependence of D is nonmonotonic; see Fig 8 .Lastly, it is found that the diffusion of the aggregates assembled with HT5 surfactants that are close to their fully extended length is smaller than that of the RMs formed with semiflexible surfactants.This is due to the high steric hindrance of the RMs made up of surfactants with high rigidity.

Conclusions
In this work, the interplay between surfactant concentration and persistence length on structural and dynamic properties of water/HT5 surfactant aggregates, forming reverse micelles as water-droplet carriers was investigated.It is found that increasing the surfactants' concentration leads to higher number of water/HT5 aggregates, although this number is generally reduced when the persistence length of the surfactants is increased.The aggregates formed are roughly spherical in shape; however, the structure of the surfactants' tail groups surrounding the RMs nuclei depends on the surfactants' rigidity.More flexible surfactants yield more compact aggregates.Our results prove that there is synergy between c s and L p , leading to enhanced diffusion of the water beads/HT5 aggregates, for certain surfactant concentrations and persistence lengths.The self-diffusion coefficient of the RMs is found to depend on the cooperative behavior between c s and L p .Also, it is noted that diffusion is improved when the number of RMs is relatively low, to ensure ample motion capacity over the entire system.Additionally, the HT5 surfactant chains must be flexible enough so that the structure enveloping the RMs nuclei is as compact as possible.Otherwise, having numerous micelles in the system leads to low RMs diffusion.On the other hand, having RMs with highly stretched external structures increases the steric hindrance of the aggregates, which in turn reduces their diffusion.These results are expected to be useful in industrial and technological applications where RMs are used.In particular, they should provide insights into processes where the use of RMs and the selection of surfactants to form such aggregates are crucial, such as in drug delivery [50], antibiotic recovery [51], protein recovery [14], and enhanced oil recovery [52].

Fig 1 (
B) shows a simplified representation of the two bonding forces (Eqs (4) and (5)) used to model the HT5 surfactant.Fig1(C) and 1(D) illustrate the consequences of having a RM made up of flexible and rigid surfactants, respectively.Eqs (1)-(5) make up the force field used here to study the

Fig 1 .
Fig 1.(a) DPD species used in all simulations, water (W, in blue), surfactant's head (H, in red), surfactant's tail (T, in green), and oil (O, in dark yellow).(b) HT5 surfactant model.The bonding and bending forces' parameters used to join consecutive beads of the surfactants are schematically represented.The constants k S and k A are the spring and angular constants of F S ij (Eq (4)) and F A ij (Eq (5)), respectively.The equilibrium distance and angle are r 0 and θ 0 ; see F S ij (Eq (4)) and F A ij (Eq (5)), respectively.In (c) and (d), one finds simplified representations of the aggregates formed by flexible and rigid surfactants, respectively.
Fig 2, comparing the effects produced by the stiffness of the HT5 on the micellization process at the lowest (c s = 6 wt%, Fig 2(A) and 2(B)) and highest surfactant concentrations (c s = 12 wt%, Fig 2(C) and 2(D)).Fig 2 shows the reduced number density profiles along the z-direction of the simulation box, omitting the nonpolar solvent profiles, for visualization purposes.All profiles displayed in Fig 2 correspond to the last simulation block of the production phase.The rest of the density profiles are found in S1 Fig in the S1 File accompanying this article.
Fig 2(A)-2(D)) and surfactants' heads profiles (red lines in Fig 2(A)-2(D)) follow each other closely.The maxima in the water beads and surfactants' heads profiles are located at the same z-position, although they differ in intensity, as a consequence of the water/surfactant molecular ratio of 2:1.Additionally, one can note that the micelle's nuclei are always covered by surfactants' tail groups, as seen in their density profiles.This is shown by the green lines in Fig 2(A)-2(D)), which always surround the maxima of the components forming the micelle's nuclei, i.e., water beads and surfactants' heads (blue and red lines in Fig 2(A)-2(D), respectively).
Fig 2(B)).This also agrees with the number of maxima of the water beads' (and surfactants' heads) density profile displayed in Fig 2(B).As c s increases, see for example Fig 2(D), it becomes increasingly difficult to associate the number of RMs formed to maxima observed in the density profile of the water

Fig 2 .
Fig 2. Density profiles of water beads (W, blue lines), surfactant's head (H, red lines), and surfactant's tail (T, green lines) beads, from the last simulation block of the systems.The left-hand panels ((a) and (c)) show the density profiles of surfactants with the lowest persistence length, while in the right-hand panels ((b) and (d)) display the density profiles of surfactants with the highest persistence length.The density profiles of the nonpolar solvent are omitted for clarity.https://doi.org/10.1371/journal.pone.0294913.g002

Fig 3 .
Fig 3. Average number of RMs, N ave RMs , as a function of the surfactant persistence length (L p ) for the four different surfactant concentrations.The averages were calculated over 18.072 ns, and the number of RMs was counted every 451.8 ps.https://doi.org/10.1371/journal.pone.0294913.g003 Fig 5 shows the RDFs between the water beads and the individual components of the surfactants, such as the head groups (W-H, Fig 5(A) and 5(B)) and the tail groups (W-T, Fig 5(C) and 5(D)).For clarity, only the RDFs for the most flexible (L p = 0.73 nm) and the most rigid surfactant models (L p = 44.99nm) are shown in Fig 5.The spatial correlation between water beads and surfactants' heads, Fig 5(A) and 5(B), show the closest spatial correlation of all RDFs, since these two components make up the core of the aggregates.Although all water beads are trapped inside the RMs nuclei, the (W-H) RDFs show closer spatial correlation than the (W-W) RDFs; see also S3 Fig in the S1 File.This is expected because the repulsive force between the water beads and the surfactants' head groups is the lowest among all the other pair of components, including those between particles of the same type; see Table1 . The information provided by the (W-T) RDFs, shown in Fig 5(C) and 5(D), confirms that the RMs nuclei are surrounded by the surfactant tail groups.This is because the (W-T) correlation is found to be displaced to larger relative distances, in comparison with the (W-H) correlation.This conclusion is supported by the vertical dashed lines in Fig 5(C) and 5 (D), which indicate that the relative position of the first maximum of the (W-T) RDFs (r = 0.94), is located at a larger relative distance than the position of the first maximum of the

Fig 4 .
Fig 4. Time evolution of the number of RMs (N RMs ) formed during the self-assembly process, as a function of time.Here, only results for the systems made up of the HT5 surfactants with the lowest (a) and highest rigidity (b) are shown, for clarity.The number of RMs was counted every 451.8 ps for a span of 18.072 ns.https://doi.org/10.1371/journal.pone.0294913.g004

Fig 5 .
Fig 5. Radial distribution functions, g(r), between water-surfactants' heads (panels (a) and (b)) and water-surfactants' tails (panels (c) and (d)) at the four different surfactant concentrations, c s .In the left-hand panels ((a) and (c)) are shown the g(r)'s of the systems made of surfactants with the lowest persistence length (L p = 0.73 nm), whereas the right-hand panels ((b) and (d)) shows the g(r)'s of the analogous systems with stiffer surfactants (L p = 44.99nm).Vertical dashed lines show the location of the first maximum of the g(r) function.https://doi.org/10.1371/journal.pone.0294913.g005 Fig 6 shows the RDFs between the surfactants' head-and tail-groups (H-T), accompanied by a snapshot of a typical RM for the most flexible (Fig 6(A) and 6(B)) and the most rigid (Fig 6(C) and 6(D)) surfactants.Fig 6 shows the influence of the surfactants' persistence length on the structure of the beads covering the RMs nuclei, revealing the average structure of the HT5 surfactants in the RMs.These results show that the arrangement of tail beads of the most flexible surfactants enveloping the RMs nuclei is more compact than those of the most rigid surfactants, whose stretched configuration is responsible for the increase in the size of the RMs.The number of maxima observed in the (H-T) RDFs, indicated by the vertical dashed lines in Fig 6(A) and 6(C), confirm the compact structure for the most flexible case.In Fig 6(A), the third maximum is barely noticeable, whereas for the most rigid case, in Fig 6(C), up to a fourth maximum can be seen, indicating a stretched surfactant configuration.

Fig 6 .
Fig 6.Averaged RDFs between surfactants' head-and tail-groups, g(r), for the most flexible ((a), L p = 0.73 nm) and the most rigid ((c), L p = 44.99nm) surfactants.Panels (b) and (d) show snapshots of typical RMs, sliced to see their core, where water contained inside the RM were removed to better visualization of the HT5 structure.The vertical dashed lines in panels (a) and (c) indicate the relative position of each maximum.R e in panels (b) and (d) refers to the average end-to-end distance of the HT5 surfactant in each case (see Eq (7) and Fig 7).https://doi.org/10.1371/journal.pone.0294913.g006

Fig 7 .
Fig 7. Mean end-to-end distance R e h i as a function of the persistence length (L p ) obtained from the average over all the surfactant concentrations.Here, the results of this work (black triangles) are compared with those for the HT9 model (blue circles) [33], and with the results of Firetto et al. [46] for the H4T4 surfactant model (red stars).The horizontal dashed lines indicate the fully extended length R FEL h i predictions for each surfactant model.The vertical and top axes are expressed in reduced units.https://doi.org/10.1371/journal.pone.0294913.g007

Fig 8 .
Fig 8. Self-diffusion coefficient (D) of the water/HT5 RMs as a function of the surfactant persistence length (L p ), for the different surfactant concentrations.The dashed lines are only guides for the eye.https://doi.org/10.1371/journal.pone.0294913.g008

Table 1 . Values of the DPD conservative force parameters, a ij (see Eq (3)), for each pair of DPD particle species.
ternary mixture systems of oil/water/surfactant, in the water/surfactant self-assembly process into RMs.

Table 3 . Values of the surfactant concentration (c s ) and the surfactant persistence length (L p ) studied in this work.
(6) values of the persistence length were obtained using Eq(6)for the different values of k A shown in this table. https://doi.org/10.1371/journal.pone.0294913.t003

Table 4 . Mean end-to-end distance ( R
e h i) of

the HT5 surfactants, for all persistence length values and surfactant concentrations modeled here.
The values of R e h i are in reduced DPD units. https://doi.org/10.1371/journal.pone.0294913.t004